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Abstract 

In  this  paper  a  method  is  proposed  for  simultaneous  measurement  of  local  groimdwater 
specific  discharge  (q)  and  Cr(VI)  mass  flux  using  a  permeable  adsorbent  device. 

The  proposed  method  allows  for  determination  of  both  q  and  as  two-dimensional 

vectors,  with  each  vector  composed  of  an  angular  (qe,  and  a  radial  (qr, 
component.  The  device  consists  of  a  cylindrical  container  whose  outer  shell  is 
constructed  of  highly  permeable  mesh  material  with  impermeable  end  caps.  The  interior 
of  the  device  is  divided  into  a  center  well  and  three  outer  sectors,  each  packed  with  a 
granular  anion  exchange  resin  having  high  adsorption  capacity  for  the  Cr(VI)  oxyanions 
Cr04^'  and  HCr04'.  The  adsorbent  in  the  center  well  of  the  device  is  loaded  with 
benzoate  as  resident  tracer.  Experiments  were  conducted  in  which  these  deviees  were 
placed  in  porous  packed  bed  columns  through  which  was  passed  a  measured  quantity  of 
simulated  groundwater  containing  Cr(VI).  The  resin  was  then  removed  from  the  four 
sectors  separately  and  extracted  to  determine  the  masses  of  Cr(VI)  and  resident  tracer 
bound  in  each.  A  linear  solute  transport  model  was  used  to  relate  the  observed  rate  of 
resident  traeer  displacement  to  qe  and  qr,  and  a  simple  mass  balance  equation  was  used  to 
calculate  Experimental  results  showed  qe  was  measured  to  an  accuracy  of  ±  15°. 
Accuracy  of  the  radial  specifie  discharge  veetor  component  depended  on  the  degree  of 
displacement  of  the  resident  tracer.  Jr*^"^  mass  flux  vector  radial  component  was  measured 
to  an  accuracy  of  ±  17%.  These  results  demonstrate  that  the  proposed  method  represents 
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a  promising  technique  for  determination  of  specific  discharge  and  contaminant  mass  flux 
in  contaminated  aquifers. 

Introduction 

Chromium  is  one  of  the  most  common  and  persistent  groundwater  pollutants  in  the 
United  States  (1).  A  considerable  amount  of  laboratory  research  has  been  devoted  to  the 
study  of  processes  that  effect  the  concentration  and  speciation  of  chromium  in  aqueous 
systems.  These  processes  include  reduction  of  hexavalent  chromium  (Cr(VI)),  the  most 
toxic  and  mobile  oxidation  state  of  chromium,  to  the  less  mobile  and  less  toxic  Cr(III)  by 
abiotic  (2,  3)  and  microbially-mediated  {4,  5,  6)  redox  processes,  adsorption  of  Cr(VI) 
oxyanions  on  metal  oxides  (7,  8,  9),  microbial  biosorption  of  Cr(VI)  {10),  oxidation  of 
Cr(III)  to  Cr(VI)  by  dissolved  oxygen  {11),  photoreduction  of  Cr(VI)  {12),  and 
precipitation  of  Cr(III)  hydroxides  {13). 

Chromium  transport  in  polluted  groundwater  has  also  been  studied  extensively 
{14, 15, 16).  However,  interpretation  of  the  results  of  these  field  studies  is  confounded 
by  the  complexity  of  Cr(VI)  chemistry  in  natural  systems  {1 7).  In  addition  to  the  above 
processes,  Cr(VI)  concentration  in  groundwater  can  be  influenced  by  dilution  due  to 
aquifer  recharge  downstream  from  a  source  zone.  Since  volumetric  recharge  rates 
usually  can  only  be  estimated  approximately,  the  effect  of  aquifer  recharge  on  Cr(VI) 
concentration  cannot  be  known  with  a  hi^  degree  of  certainty.  This  uncertainty 
prohibits  accurate  determination  of  Cr(VI)  attenuation  rates  at  affected  sites.  In  principle. 
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attenuation  rates  can  be  measured  more  accurately  using  contaminant  mass  flux  rather 
than  concentration.  The  mass  flux  J’  of  a  contaminant  i  is  the  product  of  concentration  C 
and  specific  discharge  q,  i.e.  J'  =  C'q.  (Note  that  in  this  paper,  2-dimensional  vectors  are 
indicated  by  boldface,  while  scalar  quantities  are  in  regular  typeface.)  Since  recharge 
flow  to  an  aquifer  element  increases  q  and  decreases  C  by  the  same  factor,  J'  is 
independent  of  recharge  rate.  Monitoring  Cr(VI)  attenuation  rates  through  measurement 
of  Cr(VI)  mass  flux  in  addition  to  or  instead  of  concentration  could  therefore  lead 
to  greater  accuracy  by  eliminating  the  uncertainty  in  estimated  recharge  rates.  The  ability 
to  directly  measure  could  also  be  of  significant  use  in  predicting  the  impact  of  a 
Cr(VI)  source  zone  on  a  down-gradient  drinking  water  supply  well.  The  impact  of  a 
contaminant  source  on  a  downstream  supply  well  can  be  predicted  using  a  mass  balance 
estimation  (18),  but  this  approach  relies  on  a  number  of  conditions  which  often  are  not 
satisfied.  Direct  measurement  of  and  q  across  a  transect  of  a  Cr(VI)  plume  over  a 
long  sampling  period  would  permit  more  accurate  estimation  of  future  levels  in  the 
supply  well,  by  minimizing  or  eliminating  reliance  on  assumptions  of  constant  flow  field 
and  constant  source  zone  release  rate. 

In  this  paper,  we  propose  a  general  method  for  simultaneous  measurement  of 
groundwater  q  and  J’  in  contaminated  aquifers  using  a  permeable  adsorbent  device.  In 
concept,  the  device  consists  of  a  permeable  cylindrical  shell  packed  with  a  porous  bed  of 
adsorbent  material  having  a  high  sorption  capacity  for  the  groundwater  contaminant  of 
interest.  The  ends  of  the  cylinder  are  capped  with  impermeable  material.  The  interior  of 
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the  device  is  divided  by  permeable  walls  into  four  sectors,  a  center  well  and  three  equal- 
volume  outer  sectors,  as  shown  in  Figure  1 .  A  quantity  of  a  reversibly-adsorbed  resident 
tracer  compound  is  uniformly  distributed  on  the  adsorbent  in  the  center  well,  while  the 
adsorbent  in  the  three  outer  sectors  contains  no  resident  tracer.  Sampling  is  initiated  at 
some  time  t  =  0  (Fig.  la)  when  the  device  is  inserted  into  a  sampling  well  in  the 
contaminated  aquifer.  Horizontal  infiltration  of  groundwater  through  the  device  due  to 
local  hydraulic  gradients  causes  accumulation  of  sorbed  contaminant  in  the  upgradient 
sectors  of  the  device,  while  the  passage  of  groundwater  through  the  center  well  displaces 
resident  tracer  into  the  downgradient  sectors.  After  a  sampling  period  At  (Fig.  lb)  the 
device  is  removed  from  the  sampling  well.  The  adsorbent  material  is  removed  separately 
from  each  of  the  outer  sectors  of  the  device,  extracted,  and  the  extract  fractions  analyzed 
to  determine  the  masses  of  resident  tracer  and  contaminant  adsorbed  in  each  sector  at  At. 
Groundwater  q  and  J’  can  then  be  calculated  as  two-dimensional  vectors  using 
relationships  derived  from  simple  solute  transport  models. 

The  general  method  described  above  could  be  applied  to  measurement  of  j‘^'^if  the 
appropriate  adsorbent  and  resident  tracer  were  selected.  Since  Cr(VI)  in  dilute  aqueous 
solutions  at  near-neutral  pH  is  present  predominately  as  the  oxyanions  chromate  (Cr04  ’) 
and  bichromate  (HCr04'),  and  at  higher  concentrations  as  dichromate  (CriO? ')  {19),  an 
anion  exchange  resin  would  seem  to  be  an  appropriate  choice  of  adsorbent,  and  any 
stable  anionic  compound  with  the  appropriate  partitioning  behavior  could  be  used  as  the 
resident  tracer.  The  objective  of  the  present  study  was  to  evaluate,  using  laboratory 
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(a)t  =  0  (b)t  =  At 


f  I  ff  V 

♦  ♦  + 

h  ♦  ♦  ♦ 

+  ♦  4* 

K  ♦  +  ♦ 


Adsorbent 


Adsorbent  with  bound 
resident  tracer 


Adsorbent  with  bound 
AjAjA]  groundwater  contaminant 


Figure  1.  Conceptual  diagram  showing  operation  of  a  permeable  adsorbent  device  for 
simultaneous  measurement  of  groundwater  specific  discharge  and  contaminant  mass  flux, 
(a)  Device  as  entered  in  porous  media,  (b)  Device  after  a  sampling  period  Dt. 
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experiments,  the  usefulness  of  permeable  packed  bed  devices  containing  anion  exchange 
adsorbent  for  simultaneous  measurement  of  q  and  in  porous  media.  Experiments 
were  designed  to  assess  the  accuracy  of  the  devices  in  conditions  that  simulate 
contaminated  groundwater  flow  in  an  idealized  aquifer  system.  Experimental  parameters 
were  selected  to  reveal  the  limits  of  the  device’s  useful  range  in  the  measurement  of  the 
vectors  q  and 

Experimental  Section 

Chemicals.  Benzoic  acid  (J.T.  Baker,  Phillipsburg,  N.  J.)  stock  solution  was  prepared  by 
dissolving  1.00  g  of  the  free  acid  in  10.0  mL  HPLC  grade  methanol  (Fisher  Scientific, 
Fair  Lawn,  NJ).  A  2.0  g  L'^  Cr(VI)  stock  solution  was  prepared  by  dissolving  565  mg 
potassium  dichromate  (J.T.  Baker)  in  100  mL  deionized  water.  Salt  stock  solution  was 
prepared  by  dissolving  the  following  in  1.0  L  deionized  water:  15  g  sodium  phosphate 
monobase  hydrate  (Fisher),  2.5  g  each  calcium  chloride  dihydrate  (Aldrich),  potassium 
chloride  (Sigma,  St.  Louis,  MO),  sodium  chloride  (Aldrich),  magnesium  chloride 
hexahydrate  (Fisher),  magnesium  sulfate  heptahydrate  (Sigma),  125  mg  manganese 
chloride  tetrahydrate  (Aldrich),  and  25  mg  sodium  molybdate  dihydrate  (Fisher). 
Simulated  groundwater  buffer  (SGB)  was  made  by  dissolving  5.00  g  sodium  bicarbonate 
(Aldrich)  in  about  8  L  deionized  water,  then  adding  100  mL  salts  stock  and  diluting  to 
10.0  L.  SGB  had  pH  7.0  ±  0.1,  ionic  strength  9.1  mM,  and  conductivity  0.8  mS.  Resin 
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extraction  solution  was  IM  sodium  chloride  in  deionized  water  with  pH  adjusted  to  1 1 
using  2N  NaOH.  Deionized  water  used  was  >18  MQ.  AG1-X8  anion  exchange  resin, 
20-50  mesh,  1.2  meq  ml"'  capacity,  was  obtained  in  the  hydroxide  form  from  Bio-Rad 
Laboratories  (Hercules,  CA),  and  was  converted  to  the  bicarbonate  form  by  washing 
slowly  with  12  bed  volumes  IM  NaHCOa,  then  12  bed  volumes  of  deionized  water.  The 
wet,  drained  AG1-X8(HC03  )  resin  had  bed  density  0.75  g  cm‘^,  and  the  resin 

bed  porosity  was  measured  to  be  0.50.  Benzoate  resident  tracer  compounds  were 
loaded  onto  the  resin  after  the  OH'  to  HCOb'  conversion  by  stirring  10  mL  of  resin  in  50 
mL  deionized  water,  then  adding  50  pL  2N  NaOH  and  200  pL  benzoic  acid  stock 
solution.  After  1  hour  of  stirring  the  loaded  resin  was  dump-packed  into  a  glass  column 
and  washed  with  20  bed  volumes  of  deionized  water.  Initial  adsorbed  concentration  of 
benzoate  on  the  resin  was  neither  measured  nor  carefully  controlled;  however,  the 
maximum  adsorbed  benzoate  concentration  obtainable  by  this  procedure  (assuming  all 
added  benzoate  adsorbed)  was  2  mg  mL'*,  or  0.016  meq  mL'\  less  than  2%  of  total  resin 
capacity.  Resin  portions  with  and  without  loaded  tracers  were  stored  slurried  in 
deionized  water  at  room  temperature  before  use. 

Cr(VI),  Benzoate  Isotherms.  Adsorption  isotherm  data  points  on  AG1-X8(HC03')  resin 
were  acquired  using  10  cm^  glass  vials.  To  each  tared  vial  was  added  a  portion  of  wet, 
drained  resin  weighed  to  the  nearest  mg,  and  5.0  mL  of  SGB  containing  2.0  mg  mL'*  of 
either  Cr(VI)  or  sodium  benzoate.  Duplicate  control  vials  contained  no  resin.  All  vials 
were  sealed  with  teflon-lined  septa,  placed  on  a  rotary  shaker  at  30  rpm  and  agitated  for 


Campbell,  Hatfield,  Aimable,  and  Rao  August  2001 


8 


>40h  in  a  20“C  incubator.  Aqueous  phase  concentrations  (CAq,  mg  mL’')  were 
determined  as  described  below.  Adsorbed  masses  were  calculated  from  the  difference  in 
aqueous  mass  of  solute  between  vials  with  and  without  resin.  Adsorbed  concentrations 
(Csorb,  mg  g*')  were  calculated  based  on  the  mass  of  wet,  drained  resin.  Benzoate  linear 
partitioning  coefficient  (ko^,  mL  g"')  was  calculated  from  the  slope  of  the  benzoate 
adsorption  isotherm. 

Media  Permeability.  The  hydraulie  conductivities  of  packed  beds  of  AG1-X8  (HCO3') 
resin  and  l-mm  glass  beads  and  respectively,  cm  s'*)  were  measured  using 

a  simple  constant  head  permeameter  tube  with  cross-sectional  area  Abed  =  2.41  cm  and 
length  Lbed  =  10  cm.  Steady  flow  Q  (cm^  s'*)  of  SGB  was  measured  at  three  different 
constant  head  differences  Ah  (cm)  using  a  graduated  cylinder  and  stopwatch,  then 
hydraulic  conductivity  calculated  as  K  =  QLbed/(AbedAh)  (20).  For  both  media  the  three  K 
measurements  had  relative  standard  deviations  of  <  1  %.  The  mean  of  the  three 
measurements  was  taken  as  the  hydraulic  conductivity  of  the  media. 

Device  Construction.  Permeable  adsorbent  devices  were  constructed  of  50  x  50  mesh 
milling  grade  woven  wire  cloth  (McMaster-Carr  Supply  Company,  Atlanta,  GA) 
according  to  the  design  shown  in  Figure  2.  The  wire  cloth  material  used  was  woven 
from  7.5  mil  (0.19  mm)  diameter  304  stainless  steel  wires,  had  nominal  opening  width  of 
0.3 1 8  mm,  and  was  39%  open  area.  The  outer  shell  of  the  device  was  formed  by 
overlapping  the  ends  of  an  80  x  2.6  cm  strip  of  wire  cloth.  The  overlapping  ends  were 
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Sector  dividers 


Figure  2.  (a)  Design  and  dimensions  (cm)  of  permeable  adsorbent  devices  used  in 
study,  (b)  Photograph  of  an  empty  device. 
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stitched  together  using  wire  unraveled  from  the  cloth  itself,  forming  a  cylinder  2.4  cm  in 
diameter.  The  center  well  was  made  the  same  way  using  a  4.0  x  2.1  cm  strip  of  wire 
cloth  to  form  a  cylinder  1 .2  cm  in  diameter.  The  outer  sector  dividers  were  2.1  x  0.6  cm 
strips  of  wire  cloth.  These  were  joined  to  the  center  well  at  120  degree  spacings  by 
pulling  a  wire  out  of  the  weave  along  one  of  the  long  ends  of  each  strip,  then  inserting  the 
protruding  wires  along  that  edge  into  the  openings  in  the  center  well  mesh.  End  caps 
were  cut  from  0.3  cm  -  thick  cast  acrylic  clear  sheet  (McMaster-Carr)  using  a  1-inch 
diameter  hole  saw  with  no  pilot  bit.  End  caps  cut  this  way  fit  tightly  into  the  device  outer 
shell  and  were  held  in  place  by  friction.  Devices  constructed  this  way  had  outer  diameter 
(DOev)  2  4  center  well  diameter  1 .2  cm,  and  bed  length  2. 1  cm.  Resin 
was  packed  into  the  devices  as  follows;  an  end  cap  was  pressed  into  the  outer  shell  first, 
then  the  center  well  with  attached  sector  dividers  was  inserted  into  the  shell.  The  sector 
dividers  were  arranged  into  a  radial  pattern,  then  the  outer  sectors  were  packed  with 
AG1-X8  (HCOb')  resin.  Once  the  outer  sectors  were  packed,  the  resin  served  to  hold  the 
sector  dividers  firmly  in  place.  Resin  loaded  with  benzoate  resident  tracer  was  then 
packed  into  the  center  well,  and  the  other  end  cap  was  pressed  into  place. 

Column  Experiments.  Permeable  adsorbent  devices  were  tested  in  column  systems 
configured  as  shown  in  Figure  3.  The  PVC  columns  had  inner  diameter  1 1 .5  cm 
and  length  (L^°’)  6.5  cm.  The  column  cross-sectional  area  was  104  cm^.  Columns 
were  dump-packed  with  4-mm  glass  beads  (Fisher)  to  a  depth  of  about  2  cm,  then  filled 
to  the  top  with  1-mm  glass  beads  (Supelco).  A  circle  of  woven  wire  cloth  was  used  to 
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Figure  3.  Schematic  of  porous  media  column  systems  used  to  conduct  experiments  with 
the  permeable  adsorbent  devices.  Column  dimensions  are  given  in  cm. 
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separate  the  two  bead  layers.  This  design  ensured  that  buffer  fed  through  the  bottom  of 
the  column  was  distributed  at  even  pressure  throughout  the  bottom  section  of  the  column, 
creating  an  even  upward  vertical  flow  profile  throughout  the  1-mm  bead  bed.  After 
packing,  columns  were  flushed  with  SGB  containing  5  mg  L'^  Cr(VI)  (Cr/SGB).  Flow 
through  each  column  was  siphoned  via  PVC  tubing  from  a  reservoir  which  was 
maintained  at  constant  level  by  forced  feed  pumped  from  a  main  reservoir.  Volumetric 
feed  rates  to  the  constant-level  reservoirs  were  not  precisely  controlled  and  ranged  from  3 
to  10  mL  min'^  At  t  =  0  a  permeable  adsorbent  device  prepared  as  described  above  was 
buried  in  the  glass  bead  bed  to  a  depth  of  approximately  3  cm  from  the  top  edge  of  the 
device  to  the  surface  of  the  bed.  The  device  was  placed  in  the  column  with  axis 
horizontal,  and  with  the  reference  divider  between  outer  sectors  1  and  3  rotated  an  angle 
qe  from  the  vertical.  Cr/SGB  was  then  allowed  to  flow  upward  through  the  column, 
draining  from  the  top  of  the  bed  through  weir  holes  cut  into  the  column  rim.  The  Cr/SGB 
was  collected  as  it  drained  from  the  column,  and  the  collected  volume  was  measured  in  a 
graduated  cylinder  to  the  nearest  0. 1  L.  After  a  volume  AV  (L)  had  been  collected  over  a 
sampling  period  At  (h),  flow  through  the  column  was  stopped  and  the  device  removed 
from  the  column.  The  Cr/SGB  used  in  each  run  was  analyzed  separately  to  determine 
feed  Cr(VI)  concentration  (Cf*^^  mg  L“’ ).  A  total  of  16  runs  were  conducted  in  a  4^ 
experimental  matrix,  with  the  parameters  qe  and  AV  varied  over  4  levels  each. 
Experimental  values  of  qe  were  0, 45, 90,  and  180  degrees,  while  AV  levels  were  4  -6  L, 
8  -12  L,  19  -  21  L,  and  30  -  42  L. 
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Analytical  Methods.  After  removing  the  device  from  the  column,  one  end  cap  was 
removed  from  the  device,  and  the  resin  was  scraped  out  of  each  sector  individually, 
starting  with  the  center  well,  using  a  round-end  spatula.  Care  was  taken  to  remove  as 
much  resin  as  possible  from  each  sector  without  cross-contaminating  resin  in  adjacent 
sectors.  The  resin  from  each  sector  was  transferred  to  a  50-mL  beaker,  slurried  in 
deionized  water,  then  packed  into  a  1  cm  ID  glass  column,  and  extracted  with  20  mL 
resin  extraction  solution  at  0.6  mL  min'.  This  extraction  technique  was  found  to  give 
Cr(VI)  and  benzoate  recoveries  of  100  ±  5%  (data  not  shown).  Extract  fractions  were 
assayed  for  Cr(VI)  mass  (mi^,  mg)  using  the  diphenylcarbazide  colorimetric  method 
(21).  Response  to  Cr(VI)  at  540  nm  was  calibrated  using  a  3-level  standard  curve. 
Benzoate  resident  tracer  mass  (m  mg)  in  the  extract  fractions  was  determined  by  HPLC 
using  a  25.0  x  0.46  cm  Partisil  PAC  10pm  column  (Alltech  Associates,  Inc.,  Deerfield, 
IL).  Mobile  phase  was  10%  (v/v)  HPLC-grade  acetonitrile  (Fisher),  90%  50  mM  sodium 
phosphate  pH  3.2,  at  1.25  mL  min'.  The  retention  time  of  benzoic  acid  by  this  method 
was  4.3  min.  Samples  were  diluted  a  minimum  of  1 :20  in  mobile  phase,  then  0.25  mL 
was  injected.  Detection  was  by  UV  at  220  nm.  Benzoate  UV  response  factor  was 
determined  using  a  3-level  external  standard  calibration  curve.  A  blank  and  a  calibration 
check  standard  were  analyzed  with  each  set  of  samples  to  verify  stable  response. 
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Results  and  Discussion 

Cr(VI)  and  benzoate  adsorption  isotherms  at  pH  7, 20"C  are  shown  in  Figure  4.  Cr(VI) 
partitioning  is  nonlinear,  as  is  expected  for  exchange  of  a  divalent  anion  such  as  Cr04  ' 
with  a  monovalent  anion  such  as  HCOs'  (22).  Benzoate  partitioning  is  linear  (r^  =  0.95), 
with  C®sorb  related  to  C®Aq  by  a  linear  partitioning  coefficient  (ko®,  ml  g  *) 


'^Sorb  ~  '^d'^  Aq 


[1] 


The  solid  line  in  Figure  4  represents  [1]  with  ko^  =  158  mL  g'*  for  benzoate  partitioning 
in  SGB.  Benzoate  equilibrium  partitioning  on  an  anion  exchange  resin  is  influenced  by 
the  concentrations  of  competing  anions  in  the  aqueous  phase.  Since  column  experiments 
were  conducted  with  buffer  reservoirs  open  to  the  atmosphere,  bicarbonate 
concentrations  could  have  increased  or  decreased  over  time  with  CO2  absorption  or 
release.  Therefore,  it  was  important  to  verify  that  this  ko®  value  was  valid  for  each  of  the 
16  column  experiments.  To  verify  this,  the  electrical  conductivity  of  the  Cr/SGB  feed 
buffer  was  monitored  for  each  column  experiment,  and  was  found  to  vary  by  <  8%  over 
all  runs.  Using  this  ko®  value  with  the  AG1-X8  bed  properties  and  a 
benzoate  retardation  coefficient  (R®)  was  calculated  as 

jL  B  Re  sin  [2] 

+  —  =  238 

^  Re  sin 


Hydraulic  conductivities  and  were  measured  to  be  0.39  and  0.84  cm  s’’, 
respectively.  A  hydraulic  correction  coefficient  was  calculated  as  + 
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jrBeadsy2fr Resin  ^  j  ^20)  to  accoutit  for  the  difference  between  specific  discharge  inside 

the  device  and  in  the  surrounding  media. 


Figure  4.  Cr(VI)  and  benzoate  partitioning  between  SGB  and  AG1-X8(HC03 )  resin,  pH 
7, 20'^C.  Solid  line  corresponds  to  linear  benzoate  partitioning  in  SGB  with  =  158 

mLg'^ 

Parameters  and  masses  of  Cr(VI)  and  benzoate  extracted  from  device  sectors  are 
summarized  in  Table  1  for  the  sixteen  experimental  runs.  Experimental  and  measured 
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values  of  q  and  were  calculated  as  2-dimensional  vectors  in  a  polar  (r,  0)  coordinate 
plane  perpendicular  to  the  device  axis,  with  each  vector  having  a  radial  (qr,  and  an 
angular  (qe,  component.  The  radial  vector  component  is  equivalent  to  vector 
magnitude,  while  the  angular  component  is  the  vector  direction.  For  each  run  in  Table  1, 
time-averaged  experimental  (actual)  values  of  qrAt  (m)  and  T^'^At  (g  m'^)  were  calculated 
using 


A  V 

A  Col 


[3] 


and 


yf  Af  = 


AV-C^'' 

A  col 


[4] 


Equations  to  calculate  measured  values  of  qr  and  qe  from  extracted  tracer  masses 
mi®  were  derived  from  a  simple  center  of  mass  analogy.  The  coordinates  of  the  center  of 
mass  of  a  system  of  three  point  masses  distributed  in  a  rectangular  (x,  y)  coordinate  plane 
are  given  by  Xcm  =  SmiXi/Emi,  ycm  =  Smiyi/Emi,  where  Xcm  and  ycm  are  the  coordinates  of 
the  center  of  mass,  and  Xi,  yi,  and  mi  are  the  coordinates  and  masses  of  the  three  point 
masses  (23).  These  equations  are  converted  to  polar  coordinates  using  rcm  =  (Xcm  + 
ycm^)*^  ^  and  0cm  =  tan'‘(ycm/xcm)-  We  can  treat  the  mass  of  benzoate  extracted  from  each 
sector  of  the  device  as  a  point  mass  located  at  the  geometric  center  of  the  sector.  In 
Figure  5,  the  x-  and  y-coordinates  of  the  geometric  centers  of  the  three  outer  sectors  are 
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Table  1.  Parameters  and  extracted  sector  masses  from  16  column  experiments 
using  permeable  adsorbent  devices. 

Run  Experimental  Parameters  Extracted  Sector  Masses  (mg) 

No. 


AV 

At 

qe 

. . 

Cr(Vl) 

Benzoate 

(L) 

(h) 

(deg) 

(mg  L-') 

Cr 

mi 

Cr 

m2 

rv,  Cr 
m3 

mi 

CQ 

CM 

E 

m  B 

m3 

1 

4.80 

49 

0 

4.8 

0.14 

0.44 

0.16 

0.20 

0.13 

0.20 

2 

10.3 

25 

0 

5.7 

0.30 

1.09 

0.36 

0.15 

0.05 

0.13 

3 

20.4 

76 

0 

4.8 

0.51 

1.60 

0.42 

0.15 

0.01 

0.19 

4 

30.2 

75 

0 

4.9 

0.83 

2.48 

0.96 

0.17 

0.02 

0.13 

5 

5.30 

45 

45 

4.8 

0.03 

0.39 

0.23 

0.20 

0.12 

0.20 

6 

11.7 

74 

45 

5.0 

0.06 

0.81 

0.65 

0.35 

0.16 

0.18 

7 

20.8 

115 

45 

4,8 

0.26 

1.35 

1.98 

0.41 

0.11 

0.04 

8 

41.9 

158 

45 

5.2 

0.49 

3.16 

3.33 

0.50 

0.05 

0.02 

9 

5.00 

20 

90 

4.8 

0.07 

0.24 

0.38 

0.18 

0.09 

0.12 

10 

8.30 

98 

90 

5.0 

0.20 

0.51 

0.73 

0.31 

0.24 

0.15 

11 

20.0 

95 

90 

4.8 

0.40 

0.99 

1.76 

0.29 

0.13 

0.03 

12 

40.9 

167 

90 

5.2 

0.67 

1.67 

2.51 

0.31 

0.09 

0.02 

13 

4.60 

45 

180 

4.8 

0.31 

0.32 

0.39 

0.19 

0.11 

0.08 

14 

9.80 

77 

180 

5.0 

0.75 

0.70 

0.78 

0.12 

0.25 

0.13 

15 

19.2 

114 

180 

4.8 

1.23 

0.80 

1.53 

0.07 

0.30 

0.07 

16 

39.7 

146 

180 

5.2 

1.98 

1.80 

3.50 

0.08 

0.22 

0.02 
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Figure  5.  Axial  view  of  device,  showing  rectangular  (x,  y)  coordinates  of  the  geometric 
centers  of  the  three  outer  sectors,  and  polar  reference  lines  r  =  1  and  0  =  0. 
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indicated,  along  with  the  rectangular  axes  of  the  system.  The  circle  r  =  1  is  located 
midway  between  the  eenter  well  perimeter  and  the  device  outer  wall,  and  the  line  0  =  0  is 
set  parallel  to  the  divider  between  sectors  1  and  3.  Using  this  reference  system  along 
with  the  equations  for  rcm  and  0cm  above,  we  can  now  calculate  the  polar  eoordinates  of 
the  center  of  mass  of  the  benzoate  within  the  device  (rcm^,  0cm^)  as 


Note  that  [5]  and  [6]  do  not  include  the  mass  of  benzoate  in  the  center  well.  By  locating 
the  origin  of  the  coordinate  system  at  the  device  axis,  the  coordinates  of  the  geometric 
center  of  the  eenter  well  are  (0, 0),  so  any  benzoate  mass  in  the  center  well  makes  no 
contribution  to  the  center  of  mass  calculation.  The  term  +  D^''^)/4  in  [5]  is  the 
distance  (cm)  from  the  deviee  axis  to  midway  between  the  center  well  perimeter  and  the 
outer  wall,  while  the  radical  term  caleulates  the  fraetion  of  this  distance  that  the  center  of 
mass  of  the  resident  tracer  zone  has  been  displaced.  Equation  [6]  applies  as  shown  when 
(mi®  +  m3®)/2  >  m2®;  otherwise  the  negative  of  [6]  is  used.  A  relationship  between  rcm° 
and  qr  can  be  derived  from  solute  transport  theory.  Since  benzoate  equilibrium 
partitioning  is  linear,  we  can  neglect  skewing  of  the  resident  traeer  zone  due  to  nonlinear 
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effects.  Then,  neglecting  effects  of  dispersive  transport,  the  rate  Vr®  of  benzoate 
advective  transport  within  the  device  is  related  to  qr  by 


r 


At  the  same  time,  Vr®  can  be  measured  from  rcm®  using 


yB  = 

At 


[7] 


[8] 


Equating  [7]  and  [8],  substituting  in  [5],  solving  for  the  time-averaged  quantity  qrAt,  and 
then  inserting  to  correct  for  hydraulic  conductivity  differences  gives 


J^Bf^CorrUoev  _^J^CW\  1 

0.866(wf  -/wf  ) 

2 

4- 

0.5(wf  +mf  )-ni2 

2 

4  i 

mf  +ml  +m2 

Measured  values  of  Cr(VI)  mass  flux  Jr‘^'^At  were  simply  calculated  from  extracted  Cr(VI) 
masses  using 


J5''Ar  = 


K 


Con 


Xr 


Cr 


+  WZ2 


.Cr 


Dev  T^Dev 


jUe^D 


[10] 


The  product  in  the  denominator  of  [10]  represents  the  projected  area  of  the 

device  normal  to  incident  flow. 

Figure  6  shows  qe  results  calculated  using  [6]  from  the  sixteen  sets  of  mi®  data  in 
Table  1.  Measured  q©  values  are  plotted  against  experimental  parameters.  For 
experiments  with  AV  >  8.3  L  (qrAt  >  0.83  m)  average  residual  error  of  measured  qe 


Measured  qe  (degrees) 


Campbell,  Hatfield,  Amiable,  and  Rao  August  2001 


.30  0  30  60  90  120  150  180  210 


Experimental  qe  (degrees) 


Figure  6.  Angular  component  of  specific  discharge  vector  (qe),  measured  values  from 
[6]  versus  experimental  values  for  the  16  runs.  Dashed  line  is  diagonal  for  reference. 


values  was  12  degrees,  while  the  four  experiments  with  AV  <  5.3  L  had  average  residual 
qe  error  of  50  degrees.  These  results  demonstrate  that  a  minimum  fraction  of  resident 
tracer  must  be  displaced  out  of  the  device  center  well  for  accurate  measurement  of  qe- 
Measured  qrAt  values  calculated  using  [9]  are  plotted  in  Figure  7(a)  against 
experimental  values  calculated  from  [3].  Fair  accuracy  was  obtained  for  experiments 


n/leasured  qr  (m  d"*)  Measured  qrAt  (m) 
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Figure  7.  Radial  component  of  specific  discharge  vector,  (a)  Time-averaged  specific 
discharge  (qrAt,  m),  measured  values  from  [9]  versus  experimental  values  ffom[3].  (b) 


Specific  discharge  (qr),  measured  versus  experimental  values. 
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with  AV  <  20  L  (qrAt  <  2.0  m),  with  relative  residual  error  for  those  runs  averaging  19%. 
Poorest  accuracy  was  observed  for  the  four  experiments  with  AV  >  28  L  (qrAt  >  2.8  m; 
runs  4,  8, 12,  and  16  in  Table  1,  and  circled  points  in  Figure  7),  for  which  relative 
residual  error  averaged  46%.  This  loss  of  accuracy  with  increasing  qrAt  can  be  explained 
by  using  [7]  and  [8]  to  calculate  that  for  qrAt  >  3.0  m,  rcm  >  1  -3  cm.  This  means  that  for 
those  four  experiments  a  significant  amount  of  benzoate  resident  tracer  had  begun  to 
elute  off  of  the  downstream  edge  of  the  device.  Resident  tracer  can  be  transported  within 
the  device  by  both  dispersive  and  advective  processes.  Dispersive  processes  should  have 
the  effect  of  broadening  the  resident  tracer  zone  without  displacing  its  center  of  mass. 
Note  in  Table  1  that  benzoate  was  detected  in  all  three  outer  sectors  in  every  experiment, 
including  those  sectors  that  were  upgradient  of  the  center  well.  This  observation  can  only 
be  explained  by  dispersion,  or  broadening  of  the  resident  tracer  zone  during  the 
experiments.  Dispersive  transport  of  resident  tracer  can  contribute  to  non-quantitative 
loss  of  tracer  off  of  the  device.  Non-quantitative  loss  of  resident  tracer  from  the  device 
would  increase  by  advection  as  qrAt  increases  and  by  dispersion  as  At  increases.  This 
non-quantitative  tracer  loss  causes  the  apparent  migration  of  the  tracer  zone  center  of 
mass  to  stall.  This  explains  the  observation  in  Figure  7(a)  that  measured  qrAt  values  did 
not  increase  past  3.0  m.  Figure  7(b)  is  a  plot  of  the  qr  values  obtained  by  applying  the 
trivial  calculation  qr  =  qrAt/At  to  the  qrAt  values  obtained  from  [9]  and  [3].  Excluding  the 
four  circled  points,  the  average  residual  error  in  measured  qr  values  was  19%. 
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Figure  8(a)  shows  measured  Jr‘^''At  values  caleulated  from  the  extracted  Cr(VI) 
masses  in  Table  1  using  [10]  plotted  against  experimental  Jr'^'^At  values  calculated  using 
[4].  Good  agreement  between  measured  and  experimental  values  was  observed  over  the 
sixteen  runs.  Although  Cr(VI)  was  detected  in  each  outer  sector  extract  from  every  run, 
no  Cr(VI)  was  detected  in  center  well  extract  for  any  of  the  runs  (data  not  shown).  This 
observation  shows  that  Cr(VI)  was  not  breaking  through  the  adsorbent  in  the  upgradient 
sectors.  The  good  quantitative  Cr(VI)  recovery  shown  in  Figure  8  combined  with 
observed  lack  of  breakthrough  into  the  center  well  demonstrates  that  any  Cr(VI) 
transported  into  the  device  was  effectively  bound  by  the  adsorbent.  Absence  of  Cr(VI)  in 
center  well  extracts  also  shows  that  competition  for  adsorption  sites  between  Cr(VI)  and 
benzoate  was  not  a  significant  factor  in  observed  benzoate  migration  rates.  The  masses 
of  Cr(VI)  extracted  from  downgradient  sectors  as  shown  in  Table  1  are  evidence  of 
dispersive  transport  through  the  outer  shell  of  the  device.  Figure  8(b)  shows  the  accuracy 
of  values  obtained  using  J,^''At/At  with  the  Jr‘^''At  results  from  [4]  and  [10]  for  the 

sixteen  column  experiments.  Average  residual  error  in  measured  values  was  17%. 
Using  an  equation  analogous  to  [6]  to  calculate  from  extracted  Cr(VI)  masses  gave 
poor  agreement  with  experimental  qe  values  (results  not  shown).  This  is  because  Cr(VI) 
accumulates  over  the  entire  outer  edge  of  the  adsorbent  bed,  so  that  a  small  amount  of 
quantitative  error  introduced  during  adsorbent  removal,  extraction,  and  analysis  translates 
to  large  error  in  values  calculated  using  [6].  However,  it  is  reasonable  to  assume  that 
=  qe  for  any  flow  situation  in  which  advective  Cr(VI)  transport  dominates  over  other 
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Figure  8.  Radial  component  of  Cr(VI)  mass  flux  vector,  (a)  Time-averaged  Cr(VI)  mass 
flux  (Jr^^'^At),  measured  values  from  [10]  versus  experimental  values  from  [4].  (b)  Cr(VI) 
mass  flux  measured  versus  experimental  values. 
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processes.  An  example  where  this  would  not  be  the  case  is  when  the  device  is  located  in 
a  slow  flow  field  within  a  Cr(VI)  concentration  gradient,  which  could  allow  dispersive 
transport  to  dominate  over  advection,  so  that  ^  qe-  Even  in  such  an  unusual 
situation,  however,  sampling  over  long  At  would  allow  time  for  the  gradient  to  dissipate, 
so  that  over  the  course  of  the  sampling  period  qe  would  approach  Je^’^. 

Modification  of  the  device  design  and  of  the  procedures  described  in  this  paper 
could  possibly  improve  the  accuracy  of  the  technique.  The  devices  used  in  these 
experiments  had  center  wells  one-half  the  device  outer  diameter.  This  design  gives  four 
sectors  of  equal  volume.  However,  using  a  larger  ratio  could  increase  the 

range  over  which  accurate  qr  values  are  obtained,  as  this  would  allow  the  resident  tracer 
to  migrate  farther  before  eluting  off  the  downstream  edge  of  the  device.  Using  finer- 
grained  adsorbent  and  resident  tracers  with  higher  ko  values  would  minimize  the  effects 
of  dispersive  transport  processes,  possibly  increasing  the  accuracies  of  measured  qe,  qr, 
and  values.  It  should  be  possible  to  load  a  set  of  resident  tracer  compounds  with 
various  ko  values  onto  the  center  well  adsorbent,  allowing  accurate  measurement  of  qr 
over  a  broader  range  of  qrAt  values.  An  example  of  another  anionic  compound 
appropriate  for  use  as  resident  tracer  is  pentafluorobenzoate,  which  was  found  to  partition 
linearly  in  SGB/AG1-X8(HC03')  system,  with  ko  =  401  mL/g  (data  not  shown). 

Sampling  over  long  periods  in  groundwater  with  high  qr  and  ionic  strength  would  require 
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the  use  of  a  resident  tracer  with  high  ko,  while  short  sampling  periods,  low  qr  and  low 
groundwater  ionic  strength  would  require  the  use  of  highly  mobile  resident  tracers. 

The  results  of  the  laboratory  experiments  reported  above  demonstrate  the 
potential  utility  of  permeable  adsorbent  devices  for  simultaneous  measurement  of  local 
groundwater  specific  discharge  and  Cr(VI)  mass  flux  in  contaminated  aquifers,  and  more 
generally  for  monitoring  transport  and  attenuation  rates  of  various  groundwater 
contaminants.  Further  study  is  needed  to  test  the  accuracy  and  practicality  of  these 
devices  in  actual  field  conditions. 
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